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ABSTRACT 

We investigate the destruction of dust grains by sputtering in the high- velocity 
interstellar shocks driven by supernovae (SNe) in the early universe to reveal the 
dependence of the time-scale of dust destruction on the gas density %,o in the 
interstellar medium (ISM) as well as on the progenitor mass M pr and explosion 
energy E 5 i of SN. The sputtering yields for the combinations of dust and ion 
species of interest to us are evaluated by applying the so-called universal relation 
with a slight modification. The dynamics of dust grains and their destruction by 
sputtering in shock are calculated by taking into account the size distribution of 
each dust species, together with the time evolution of temperature and density 
of gas in spherically symmetric shocks. The results of calculations show that the 
efficiency of dust destruction depends not only on the sputtering yield but also on 
the initial size distribution of each grain species. The efficiency of dust destruction 
increases with increasing E 5 i and/or increasing n^ , but is almost independent 
of M pr as long as E§i is the same. The mass of gas swept up by shock is the 
increasing function of E$i and the decreasing function of n^ Q . Combining these 
results, we present the approximation formula for the time-scale of destruction 
for each grain species in the early universe as a function of £"51 and tih,o- This 
formula is applicable for investigating the evolution of dust grains at the early 
epoch of the universe with the metallicity of Z < Z & . The effects of the 
cooling processes of gas on the destruction of dust are briefly discussed. 
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1. Introduction 

Population III (Pop III) stars formed from metal-free gas clouds where H 2 molecule is 
the main coolant of gas are considered to be more massive than 10 2 M (Bromm & Larson 
2004 and references therein), though there are still numerous uncertainties regarding the 
mass range of the first stars due to a poor understanding of the relevant radiative feedback 
in the late accretion phase (see for example Tan & McKee 2004). Once the primordial gas is 
enriched by dust grains, dust have great influences on the subsequent formation and evolution 
history of stars and galaxies; the star formation rate (SFR) in the metal-poor star-forming 
clouds is enhanced via the formation of H2 molecules on the surface of dust grain (Hirashita 
& Ferrara 2002; Cazaux & Spaans 2004). In addition, the cooling of gas by thermal radiation 
from dust themselves makes even the gas clouds with metallicity as low as 10~ 5 Z & fragment 
into low-mass gas clumps (Schneider et al. 2002, 2003; Omukai et al. 2005). Thus, dust 
grains in the early universe may cause the characteristic mass scale of stars to shift from very 
high mass to low mass of ~ 1 M & observed at present time and strongly affect the evolution 
of the initial mass function (IMF) at the early epoch of the universe. 

Also, dust grains in the early universe play crucial roles in revealing the structure and 
evolution of the universe from the relevant observations because they absorb stellar light and 
reemit it at infrared (IR)-submillimeter (submm) wavelengths. The thermal emission from 
dust at high redshifts can distort the cosmic microwave background (CMB) radiation (Loeb 
& Haiman 1997), and the obscuration and reddening of starlight by dust residing in the 
early interstellar and/or intergalactic space lead to the serious underestimate in evaluating 
the cosmic SFR from the observations directed toward the higher redshifts (e.g., Hauser & 
Dwek 2001). 

In fact, the recent observations have confirmed the presence of large amounts of dust 
grains at high redshifts of z > 5. The submm (Priddey et al. 2003; Robson et al. 2004) and 
millimeter (Bertoldi et al. 2003) observations reported continuum thermal emission of dust 
grains from some quasars at z > 5, and have suggested a large amount of cold dust reaching 
up to 10 8 -10 9 M in the host galaxies of quasars. Maiolino et al. (2004a) discovered the 
substantial extinction of stellar light by dust in the broad absorption line (BAL) quasars 
at z > 4.9. They have shown that their extinction curves could be different from those of 
low-redshift quasars at z < 4, and have suggested that the difference may reflect different 
formation and evolution mechanisms of dust grains at z > 5. 
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The main sources of dust grains in the early universe at z > 5 are considered to be 
supernovae (SNe) evolving from the massive stars because of their short lifetime. Dust 
grains formed in the ejecta of SNe are injected into the interstellar medium (ISM), and are 
subjected to destructive processes in the interstellar shocks induced by the ambient SNe 
(e.g., Jones et al. 1994; Dwek, Foster, & Vancura 1996). Thus, the size distribution and 
amount of dust in the early interstellar space are determined by the balance between their 
production in SNe and their destruction in shocks, and change with time according to the 
star formation activity (e.g., Dwek 1998). 

How much amount of dust grains are destructed in shocks as well as how much dust 
grains absorb stellar light and reemit it by thermal emission heavily depends on their chemical 
composition, size and amount. Thus, it is essential to clarify the evolution of dust in the early 
universe by treating the formation and destruction processes of dust grains in a consistent 
manner, in order to elucidate the SFR and the IMF during the early evolution of the universe 
from the future observations. Nozawa et al. (2003) investigated the chemical composition, 
size distribution and amount of dust grains formed in the ejecta of Pop III SNe. In this 
paper, as the second step to reveal the evolution of dust grains in the early universe, we 
explore dust destruction in the high-velocity interstellar shock driven by the SN explosion, 
adopting the dust model by Nozawa et al. (2003) for the initial dust residing in the early 
ISM. 

We focus on dust destruction in nonradiative shock where dust grains are considered to 
be efficiently destructed by non-thermal and thermal sputtering (e.g., McKee 1989) because 
the nonradiative shock has a high shock velocity (> 100 km s _1 ) and a high gas temperature 
(> 10 6 K). However, only a few set of experimental data on sputtering yield are available 
for the projectile-target combinations of astrophysical interest. Hence, we first evaluate the 
sputtering yields for the combinations of dust species produced in Pop III SNe and ion species 
of interest to us, collecting a large amount of sputtering data and applying the universal 
relation for sputtering yield proposed by Bohdansky (1984) with a slight modification so 
as to reproduce well the experimental data. Then, we investigate the dependence of the 
efficiency of dust destruction for each grain species not only on the gas density nH,o i n 
the ISM but also on the explosion energy E 5 i and progenitor mass M pr of SN driving the 
interstellar shocks by applying the hydrodynamical models of Pop III SNe by Umeda & 
Nomoto (2002). Finally, we evaluate the time-scale of destruction for each grain species. 

The efficiency of dust destruction by sputtering depends on the temperature and density 
of gas and the velocity of dust relative to gas in shock as well as the sputtering yield and 
the dust size. Thus, in the calculation, we solve the time evolution of temperature and 
density of gas in spherically symmetric shocks, including the cooling of gas by the atomic 
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process, the inverse Compton scattering and the thermal emission of dust. In coupled with 
the hydrodynamical calculations for gas, we carefully treat the erosion of dust by thermal 
and non-thermal sputtering caused by motion of dust relative to gas, taking into account 
the size distribution of each dust species. 

This paper is organized as follows: We first define the efficiency of dust destruction in 
§2. In § 3, we describe the properties of the ISM in the early universe together with the 
dust model adopted in the calculations, and in § 4, we present the basic equations of the 
hydrodynamic calculations for gas and the method of the simulation for dust destruction 
in shock. The sputtering yield of each grain species necessary for the calculations of dust 
destruction is evaluated in § 5, and the grain physics in shock is described in § 6. The results 
of calculations of dust destruction are given in details for the model of interstellar shock 
that is driven by a SN with E 51 = 1 and M pr = 20 M and is propagating into the ISM 
with hh,o = 1 cm -3 , and the approximation formula for the time-scale of dust destruction 
as a function of E 5i and nn,o is derived in § 7, where the effects of the cooling of gas on 
the evolution of nonradiative shock and the destruction of dust are briefly discussed. The 
summary is presented in § 8. 



2. The definition of the efficiency of dust destruction 

Our main goal is to reveal the time-scale of dust destruction by interstellar shocks in 
the early universe. In general, the time-scale of destruction tsnj for grain species j is defined 
as (e.g., McKee 1989) 



— 1 -^swept . , 



where ej is the efficiency of destruction of grain species j, M swept is the mass of gas swept 
up by shock until the shock velocity V shock decelerates below 100 km s -1 , M ISM is the mass 
of gas and dust in the ISM, and 7 S n is the effective SN rate. McKee (1989) has estimated €j 
to be nearly constant and considered that M swept is proportional to the SN explosion energy 
£51 in units of 10 51 erg. However, it is expected that both ej and M swept depend on not 
only £51 but also the progenitor mass M pr of SN and the number density of gas nn,o in the 
ISM. Hence, we aim at deriving the dependences of tj and M swcpt on E51, M pr and nn,o- The 
truncation time t tr being defined as the time at which V s hock drops down to 100 km s _1 , the 
efficiency of destruction tj for each dust species is calculated by 

= Mjf(t = t tI ) 

3 M^\t = t tr y [) 
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where M^ st (t) and M^ ept (t) are, respectively, the mass of dust destructed and the mass 
of dust swept up by shock, and the mass of gas M swcpt swept up by spherically symmetric 
shock is given by 

An 

M swcpt ~ yPo^shock^ = *tr), (3) 
with the gas density p in the ISM and the travel distance of the shock front R s h oc k(t)- 



3. Properties of the ISM in the early universe 

3.1. Density, temperature, and metallicity of the ISM 

We assume that the ISM in the early universe is homogeneous and stationary (v — 0), 
referring to the results of the simulations on the formation of the early H II region by the 
first stars (Kitayama et al. 2004), which have shown that the radiative feedback from the 
massive stars causes the ambient ISM to be homogenized with the mean gas density of < 
1 cm~ 3 . In the calculations, we consider the hydrogen number density in the ISM to be 
in the range of 0.03 cm -3 < nn,o < 10 cm -3 to reveal the dependence of the efficiency of 
dust destruction on the gas density in the ISM. The ISM surrounding the SN explosion is 
considered to be ionized by the radiation of the massive progenitor stars and heated to by 
T ~ 10 4 K (Kitayama et al. 2004; Machida et al. 2005), and then T is expected to decrease 
below 10 4 K where the recombination rate can be high even for the metallicity as low as 
Z = 10~ 4 Z Q (Wolfire et al. 1995). However, the evolution of gas temperature in postshock 
flow is almost independent of To until the end of the Sedov-Taylor phase of high-velocity 
shock considered in this paper, and the temperature of gas T s h oc k at the shock front does 
not decrease below several of 10 4 K at the truncation time, as is shown in the result of 
calculations (see 7.1). Thus, we take T = 10 4 K as the representative gas temperature in 
the ISM, independently of the density of gas. 

The truncation time as well as the time evolution of gas temperature and density struc- 
tures in postshock flow is considered to be affected by the cooling of gas depending on the 
metallicity of the ISM. Especially, when T shock < 10 5 K, the efficient line coolings of H and 
He ions cause the shock velocity to decelerate rapidly, and resultingly regulate the trunca- 
tion time. However, the gas cooling function is almost independent of the metallicity for 
Z ^ -Z© corresponding to [Fe/H] < —3 since the metal lines do not contribute to the 
cooling of gas at T > 10 4 K (Sutherland & Dopita 1993), and we adopt the gas cooling 
function for the zero metal case given in Sutherland & Dopita (1993) (see 4.1). Note that 
the efficiency of dust destruction, which is defined by the mass ratio of dust destructed to 
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dust swept up by shock, does not depend on the metallicity as long as Z < 10~ 3 Z Q , despite 
the fact that the amount of dust in the ISM is proportional to the metallicity, and in the 
calculations we assume that the metallicity of the ISM in the early universe is Z = 10~ 4 
Z & to investigate the destruction of dust in the ISM. The elemental composition of gas in 
the ISM is given in Table 1 where the C, N, and O metals are included so as to adjust the 
metallicity to Z = 10~ 4 Z Q and their abundances are simply evaluated by multiplying the 
metallicity and the primordial ratios tabulated in Sutherland & Dopita (1993); note that 
these metals play no role in the physical processes relevant to dust destruction because of 
the much lower abundances than H and He. The electron number density n e is fixed to 1.128 
n H given by Sutherland & Dopita (1993). 



3.2. Model of dust in the early universe 

Dust grains in the ISM at the early epoch of the universe considered in this paper 
are expected to be predominantly produced in the ejecta of SNe. To reveal what kind of 
dust grains form in the early universe, Todini & Ferrara (2001) calculated dust formation in 
primordial core-collapse supernovae (CCSNe), assuming the uniform elemental composition 
and gas density within the He core. They have shown that the newly formed grains are 
amorphous carbon with size around 300 A, and Al 2 03, MgSiC>3, Mg 2 Si04 and Fe3C>4 whose 
sizes are around 10-20 A. By extending the model of dust formation by Todini & Ferrara 
(2001) to pair-instability supernovae (PISNe) whose progenitors evolve from metal-free stars 
with M pr = 140-260 M (Umeda & Nomoto 2002; Heger & Woosley 2002), Schneider, Ferrara, 
& Salvattera (2004) found that the typical sizes of dust formed in PISNe range from 0.001 
to 0.3 /im depending on grain species. The dust models by Todini & Ferrara (2001) can 
successfully reproduce the extinction curve observed toward a BAL quasar SDSS1048+46 
at z — 6.2 (Maiolino et al. 2004b) and the IR spectral energy distribution (SED) of blue 
compact dwarf galaxy SBS 0335-052 (Takeuchi et al. 2003) except for the slight overestimate 
of the far IR continuum; SBS 0335-052 is a good target in the local universe suitable for 
scrutinizing the properties of dust formed in SNe at high redshifts because of its young age 
(< 10 7 yr) and low metallicity (1/41 Z Q ) (Vanzi et al. 2000). 

However, as discussed by Kozasa, Hasegawa, & Nomoto (1989), the species of dust 
formed in the ejecta largely depend on the elemental composition in the He core, and their 
sizes are affected by the time evolution of gas density and temperature. Thus, Nozawa 
et al. (2003) extensively investigated dust production in Pop III CCSNe and PISNe, by 
considering two extreme cases for the elemental composition in the ejecta; the unmixed case 
with the original onion-like structure and the uniformly mixed case within the He core. In 



-7- 



the calculations of dust formation, the radial profile of gas density is properly treated and 
the time evolution of gas temperature is calculated by solving the radiative transfer equation 
including the energy deposition of radioactive elements, which are not taken into account by 
Todini & Ferrara (2001). Nozawa et al. (2003) found that a variety of grain species (C, Si, 
Fe, FeS, AI2O3, MgSiC>3, Mg 2 SiO,4, Si0 2 and MgO) condense in the unmixed ejecta reflecting 
the difference of the elemental composition at the formation site, and that only oxide and 
silicate grains (A1 2 3 , MgSi0 3 , Mg 2 Si0 4 , Si0 2 and Fe 3 4 ) form in the mixed ejecta with 
C/O < 1. The average sizes of these grains span the range of 0.001-1 /im depending on 
the elemental composition and gas density at the location of formation in the ejecta. The 
species of dust formed in Pop III SNe and the behaviors of their size distributions are almost 
independent of the progenitor mass, but the relative mass fraction of each grain species is 
dependent on. 

Based on the dust models by Nozawa et al. (2003), Hirashita et al. (2005) have shown 
that the extinction curve of a quasar SDSS1048+46 at z — 6.2 can be better explained by 
the dust grains produced in the unmixed Pop III CCSNe than those in the mixed Pop III 
CCSNe. Takeuchi et al. (2005) have also shown by using the dust models for M pr = 20 M© 
that only dust grains formed in the unmixed ejecta but not in the mixed ejecta can give an 
excellent fit to the observed SED of SBS 0335-052 over near to far IR wavelengths. These 
works strongly suggest that the comparisons with the observations relevant to dust in the 
early universe prefer dust grains produced in the unmixed Pop III CCSNe rather than those 
in the mixed Pop III CCSNe. 

Therefore, as the representative model of the initial dust grains residing in the early 
ISM, we adopt the size distribution and abundance of each grain species obtained by the 
calculations of dust formation in the unmixed ejecta of Pop III SN with M pr = 20 M by 
Nozawa et al. (2003). For comparison, we also consider the model of dust formed in the 
mixed ejecta of the same SN. Hereafter, the models of dust formed in the unmixed and mixed 
ejecta are referred to as the unmixed grain model and the mixed grain model, respectively. 
The initial values of the depletion factor defined by the ratio of the total dust mass to the 
total metal mass produced in the SN is 0.224 for the unmixed grain model and 0.297 for the 
mixed grain model. The bulk density p^j, mass fraction Aj and average size a ave j of each 
grain species are summarized in Table 2, and their size spectrums are given in Figure 6a for 
the unmixed grain model and in Figure 7a for the mixed grain model, where the bin width 
is 0.1 dex. 

It is useful for the following discussion to address here the behavior of size distribution of 
each dust species produced in the unmixed ejecta. C, Fe and Si0 2 grains have a lognormal- 
like size distribution with the relatively large average radius of > 0.01 /im, while the size 
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distribution functions of FeS, Mg 2 Si0 4 and MgO grains are approximately power-law. The 
average radii of A1 2 3 and MgSi0 3 grains are very small (< 5 x 10~ 3 /im). Note that Si 
grain shows a bimodal size distribution; the grains with the size less than 0.03 occupy 
almost 100 % in number but only 0.2 % in mass. Thus, as the average radius of Si grain, 
we take 0.25 /im calculated for the grains larger than 0.03 /xm including 99.8 % in mass. 
The size distribution function summed up over all grain species is well fitted with a broken 
power-law formula whose index is —3.5 for the radius larger than 0.06 /im and —2.5 for the 
smaller one. 



4. Physics of interstellar shock 

The deceleration rate due to drag force and the erosion rate due to sputtering of dust 
grains depend on the time evolution of gas temperature and density in postshock flow. In 
this section, we present the basic equations and initial conditions for the hydrodynamical cal- 
culation of the time evolution of shock and the method of the simulation for dust destruction 
in shock. 



4.1. The basic equations for the gas dynamics 

The hydrodynamic equations for gas to calculate the time evolution of spherically sym- 
metric interstellar shock driven by a SN explosion, using conventional physical variables, are 
given by 

t + ^VH0 (4) 

d Id dP 

clTJ 1 8 

+ - 2 Q- r [r 2 (U + P)v]=-A(n H ,n (6) 

where P = pkT j [im^ is the gas pressure with the mean molecular weight /i and the hydrogen 
mass ran, and the sum of kinetic energy and thermal energy per unit volume is 

with adiabatic index 7, and we adopt 7 = 5/3 in the calculation. 
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The cooling rate of gas A(nn, T) in units of erg s 1 cm 3 is expressed by 
A(n H , T) = n e n H A gas (T) + A ic (T) + A d (n H , T), 



8 



where each term represents the rate of cooling by the atomic process, by the inverse Compton 
scattering, and by the thermal emission of dust. In the calculations, we adopt the gas cooling 
function A gas (T) for the zero metal case given in Sutherland & Dopita (1993) and ignore the 
contribution of cooling by metal ions ejected from dust by sputtering, because as mentioned 
in 3.1, the metals do not contribute to the cooling of gas for the metallicity of Z = 1CT 4 
Z Q (Sutherland & Dopita 1993). For the inverse Compton scattering caused by collisions 
of hot electrons in shock with CMB photons, the cooling rate A ic (T) is given by (Ikeuchi & 
Ostriker 1986) 



where T 4 is the gas temperature in units of 10 4 K, and we take z = 20 in the calculation. 
The cooling function A d (nn, T) through the thermal emission from dust heated by collisions 
with gas particles in shock is calculated by (Dwek 1987) 



where n^j and fj(a) are, respectively, the number density and size distribution function of 
dust species j at a given t and r. The collisional heating rate Hj(a,T,nu) of grains species 
j is given in 6.3. 

The time evolution of gas temperature and density in shock driven by SN is characterized 
by the explosion energy E 5 i and the progenitor mass M pr which specify the structures of 
density and velocity of gas in the ejecta. Thus, as the initial condition for interstellar shocks, 
we adopt the hydrodynamical models of the ejecta of Pop III SNe with various E 51 and M pr by 
Umeda & Nomoto (2002) to investigate the dependence of the efficiency of dust destruction 
on E 5 i and M pr . Note that in the early universe with Z = 10~ 4 Z & , the progenitors of SNe 
are not only Pop III stars but also massive stars of the subsequent generation. However, 
the structures of density and velocity in the SN with a given E 51 and M pr do not so much 
depend on the initial metallicity for the progenitor star whose metallicity is less than 10~ 4 
Zq (Umeda, H. 2002 in private communication). In the calculations, we consider ordinary 
CCSNe with the explosion energy of E^i—l, hypernovae (HNe) whose explosion energy is 
more than ten times that of CCSNe, and PISNe whose progenitor is very massive (M pr = 150, 
170 and 200 M ), and hereafter CCSNe and HNe with progenitor mass of M pr = 13-30 M© 
are referred to as Type II SNe (SNe II). The ejecta models used in the calculations are given 
in Table 3, where the labels C, H and P represent CCSNe, HNe and PISNe, respectively, 
and the numerical value denotes the mass of the progenitor in units of solar mass. 



A ic (T) = 5.41 x 10" 32 (1 + z) 4 n e T A 



(9) 




(10) 
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4.2. The method of simulations 

Equations (4)-(6) are solved by using the flux-splitting method with second-order accu- 
racy in space and first-order accuracy in time (van Albada, van Leer, & Roberts 1982; Mair 
et al. 1988). This algorithm is one of the upwind schemes for the Euler equations and is 
well suited to solve problems involving a shock. 

The calculations for the time evolution of shock and dust destruction are performed 
via the following two steps. In the first step, we assume that the SN ejecta collides with 
the ISM in ten years after the explosion and calculate the time evolution of gas density and 
temperature without including dust, based on the ejecta models of Umeda & Nomoto (2002). 
The spatial mesh number is 500-5000 depending on the gas density in the ISM, and the inner 
20 meshes are assigned to the ejecta. Note that dust destruction in this first step does not 
contributes to the evaluation of the destruction efficiency because the mass of gas swept up 
by shock is less than 0.1 % of that by t = t tT . Also, the cooling of gas by dust has no effects 
on the subsequent structure of shock. In the second step, fixing the spatial mesh number 
to 430, we rearrange the density, temperature and velocity of postshock gas obtained from 
the first step into the inner 30 linear meshes. Then, including dust, we calculate the time 
evolution of shock as well as dynamics and destruction of dust. 

By treating dust as a test particle, the physical processes relevant to dust are calculated 
as follows. Consider a stationary (vd = 0) dust grain characterized by its composition 
and number density with radius between a and a + da at a given position Rd in the ISM. 
Once the grain enters into the shock front, the relative velocity to gas Wd = v — Vd is 
assigned. Then, for the temperature T and density n H of gas at the shock front, we calculate 
the deceleration rate dw^/dt due to drag force, the erosion rate da/dt by sputtering and 
the heating rate H(a,T,nu) by collisions with gas. By updating the relative velocity by 
w d = Wd + (dwd/dt)At and the grain radius by a' = a + (da/dt) At for a given small time 
step At, the velocity and position of the dust at time t' — t+ At are calculated by v ' d = v+w' d 
and R' d = Rd + v' d At. The procedure described above is repeated for the gas temperature 
T and density rin at the position R' d of the grain. Note that the dust grain is considered to 
be completely destroyed if the grain radius is smaller than the nominal monomer radius of 
condensate given in Nozawa et al. (2003). 

5. Evaluation of sputtering yield 

Sputtering is divided into two categories reflecting the difference in the underlying mech- 
anisms; physical sputtering and chemical sputtering. Physical sputtering invokes a transfer 
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of kinetic energy from the incident particle to target atoms and subsequent ejection of those 
atoms which have acquired kinetic energy enough to overcome the binding forces exerted 
by the target. Chemical sputtering invokes a chemical reaction induced by the impinging 
particles which produces an unstable compound at the target surface (Sigmund 1981). In 
the high-velocity shocks considered in this paper, dust grains are predominantly destroyed 
by physical sputtering because dust temperature is at most 100 K in the postshock region 
(Dwek 1987); chemical sputtering is of no importance since this process is realized at the 
temperature as high as 800 K (Roth 1983). Thus, in what follows, we simply refer to physical 
sputtering as "sputtering". The efficiency of destruction by sputtering is characterized by 
the sputtering yield defined as the mean number of emitted atoms per incident particle. The 
sputtering yield depends on the impact energy of a projectile and the incident angle to the 
target surface as well as on the projectile-target combination. However, only a few set of 
experimental data on sputtering yield are available for the combinations of dust materials 
and ion species of interest to us. 

One method evaluating the sputtering yields for the projectile-target combinations of 
astrophysical interest is to apply the Monte-Carlo code for the TRANSPORT OF IONS IN 
MATTER (TRIM, Ziegler, Biersack, & Littmark 1985) that simulates the energy loss of a 
projectile ion with given impact energy and incident angle in a solid (Field et al. 1997; May 
et al. 2000; Bianchi & Ferrara 2005). The TRIM code requires as input parameters not 
only the surface binding energy of the target, but also the displacement energy and the 
lattice binding energy E h . Particularly, uncertainties in the value of E d are a significant 
source of error in the calculated sputtering yield (May et al. 2000; Bianchi & Ferrara 
2005). In fact, the sputtering yields computed by this simulation code have been known to 
significantly differ from the experimental data at low impact energies (Bianchi & Ferrara 
2005). Furthermore, it is time-consuming to combine the Monte-Carlo simulation with the 
hydrodynamical calculation for dust destruction. 

Another one is to apply the universal relation which is the analytic formula derived by 
Bohdansky (1984) for the energy dependence of sputtering yield at normal incidence for a 
monatomic solid. To investigate the destruction of interstellar dust such as C, Fe, SiC, Si02 
and H 2 ice, Tielens et al. (1994) have applied the universal relation by treating a constant 
K in this formula as a parameter, and determined it by comparing the available experimental 
data with the sputtering yield calculated by the universal relation. They have shown that 
the universal relation can be applied to not only monatomic solids but also compound solids 
by taking the appropriate average values for the atomic and mass numbers and the surface 
binding energy under the assumption that each of element consisting of the target is sputtered 
off with the same rate. However, the agreement of the calculated sputtering yield with the 
experimental data is not so much good for the ion-target combinations with the mass ratio 
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of target atom to incident ion ranging from 0.7 and 2. 

Therefore, we apply the universal relation for sputtering yield with a slight modification 
so as to reproduce well the experimental data, collecting and examining minutely a large 
amount of sputtering data, in order to evaluate the sputtering yields for the combinations of 
dust and ion species of interest to us. Following Tielens et al. (1994), we treat K as a free 
parameter and determine it by fitting to the experimental data for the grain species with the 
yield data covering a wide range of impact energy of the projectiles. For the grain species 
for which little or no sputtering data are available, we deduce the value of K by fitting to 
the sputtering yield calculated by means of the Monte-Carlo code for the EROSION AND 
DEPOSITION BASED ON DYNAMIC MODEL (EDDY, Hirooka 2001). The EDDY code 
requires only the surface binding energy as an input parameter and allows us to predict the 
sputtering yield for specific impact energy and incident angle more easily than the TRIM 
code, though the energy range applicable to the calculations is limited to 0.1-10 keV. In 
the following subsection, we at first introduce the universal relation for sputtering yield, and 
then determine the value of K for each grain species by fitting to the available experimental 
data and/or the results of the EDDY simulation. 



5.1. Universal relation for sputtering yield 



For backward sputtering and with the neglect of inelastic energy losses, the sputtering 
yield at normal incidence Y®(E) by the projectile i impacting with energy E is given by 
(Bohdansky 1984) 



Y°(E) = 4.2 x 10 
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where U$ is the surface binding energy in units of eV, //j = Mj/Mj is the ratio of the mass 
number of the target atom to the incident ion, ctj is the energy-independent function of Hi, 
and K is considered to be a free parameter to be adjusted to reproduce the experimental 
data (Tielens et al. 1994). The threshold energies E t h have been obtained by fitting the 
yield data for low-energy sputtering and are approximately given by (Bohdansky et al. 1980; 
Andersen & Bay 1981) 
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where gi = AMiM d /(Mi + M d ) 2 is the maximum fractional energy transfer possible in a 
head-on elastic collision. The function Si(E) is the nuclear stopping cross section in units of 
ergs cm 2 and can be expressed by the following universal relation (Sigmund 1981) 

Si(E) = 4na sc ZiZ d e 2 M ^ l M ^ s i (e i ), (13) 

where e is the elementary charge, and Zi and Z d are the atomic numbers of the projectile 
and the target, respectively. The screening length a sc for the interaction potential between 
the nuclei is given by 

a sc = 0.885a [z] /3 + Z d 2/3 ) (14) 

with the Bohr radius a = 0.529 A. The function Sj(ej) depends on the detailed form adopted 
for the screened Coulomb interaction and can be approximated by (Matsunami et al. 1980) 

«<«) = 3.441^ + 2.718) (15) 

1 + 6.35 y/Ti + a (-1.708 + 6.882^] 



where the reduced energy e« is given by 

Md ° sc ,E. (16) 



M t + M d Z t Z d e 2 ' 

The value of a.i{ni) in Equation (11) depends on how the distribution of energy deposited 
in the target is approximated. Sigmund (1969) suggested the three types of approximations 
according to the adopted distribution function (Gaussian, corrected Gaussian, and non- 
Gaussian). The values of a^A**) obtained by these approximations are very similar (Sigmund 
1969) and is approximately expressed for fa < 10 as follows (e.g., Bohdansky 1984); 

ai = 0.2 for fa < 0.5 

Oi = 0.3fa /3 for 0.5 < fa < 10. (17) 

However, Equation (17) for ctj slightly overestimates the sputtering yields for the projectile- 
target combinations with the mass ratio of 0.7 < fa < 2. Hence, from a comparison with a 
large amount of sputtering data, we apply the following formula for c^; 

an = 0.2 for fa < 0.5 



a 



= 0.1 fa 1 + 0.25 (fa - 0.5) 2 for 0.5 < fa < 1 

= 0.3 (fa - 0.6) 2/3 for 1 < fa. (18) 
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As we shall show below, this modified function ai can give an excellent agreement with the 
sputtering data for any values of fa considered here (0.3 < fa < 56). 

By being K as a parameter, the surface binding energy U must be specified to evaluate 
the normal-incidence sputtering yield Y^{E) from Equation (11). We assume U to be equal 
to the sublimation energy and evaluate it from the JANAF Thermochemical Tables (Chase 
et al. 1985) except for C and AI2O3 whose values of Uq are discussed later. The values of Uq 
for each grain species are summarized in Table 4 along with the values of Z d and M d . 

5.2. Sputtering yield of each grain species 

In order to evaluate the sputtering yield for the projectile-target combinations of interest 
to us, we determine the value of K by fitting to the experimental data for C, Si, Fe, Si0 2 
and AI2O3 grains for which a sufficient amount of sputtering data are available. For MgO, 
FeS and Fe3C>4 grains for which little or no yield data are available, we evaluate the value 
of K by fitting to the results of the EDDY simulation. The experimental and simulated 
sputtering data at normal incidence, for which the references are summarized in Table 4, 
are given in Figure 1 along with the best-fitting theoretical yield calculated by the universal 
relation (solid lines). In Figures la-le, the experimental data are shown for the following 
ion species; H + (open circles), D + (asterisks), He + (open squares), Ne + (crosses) and Ar + 
(open triangles). Even in the case that the experimental data can be sufficiently available, 
for comparison, we also plot the results of the EDDY calculations for the projectiles H + 
(filled circles), He + (filled squares) and Ar + (filled triangles). 

As for C, the sublimation energy estimated to be Uq = 7 A3 eV from the thermodynam- 
ical data cannot reproduce the experimental data. Tielens et al. (1994) have shown that 
the ion fluences efficiently amorphize the surface layer of carbon material and reduce the 
surface binding energy. Thus, we take Uq = 4.0 eV following Tielens et al. (1994), and this 
value gives a good fit to the measured yields with K = 0.61 (Fig. la) being different from 
K = 0.65 by Tielens et al. (1994); the difference stems from the adopted formula for ctj. 
The sputtering yield calculated by the EDDY code is enhanced for H + by a couple of factor 
compared to the experimental data, but gives a good agreement for He + . 

The sputtering data for Si and Fe targets are sufficiently available. The modified uni- 
versal relation employing Equation (18) for ctj shows the excellent agreement with the data 
points for all projectiles considered here (Figs, lb and lc), by adopting K = 0.43 and 
K = 0.23, respectively. The EDDY calculations for Si and Fe targets also present a good 
agreement with their measured data. The dotted lines indicate the sputtering yields by 
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Ar + projectile calculated by using the values of K derived above but adopting aii given 
by Equation (17). It can be seen from these figures that Equation (17) overestimates the 
sputtering yields for the projectile-target combinations with the mass ratio of 0.7 < /Xj < 2 
as mentioned in 5.1. It should be emphasized that the previous formula cannot reproduce 
the experimental data better than the modified universal relation whatever values of K are 
selected. 

A very good fit of the universal relation to the experimental data is obtained for Si0 2 
target by taking K = 0.1 (Fig. Id). This justifies the application of the universal relation to 
compound materials as shown by Tielens et al. (1994). The results by the EDDY simulations 
are consistent with the experimental data, but a comparison with the theoretical curves of 
sputtering yield shows a little enhancement for H + and Ar + projectile at E < 200 eV. 
For MgSi0 3 and Mg 2 Si0 4 materials for which the sputtering data cannot be found in the 
literatures, as is the same as Tielens et al. (1994), we assume K to be 0.1, considering 
Si0 2 as a representative of silicates. In fact, this could be reasonable since K ~ 0.1 for the 
compounds including oxygen atoms (see Table 4). 

Given the surface binding energy of U = 6.37 eV evaluated from the thermodynamical 
data, the yield data of AI2O3 cannot be reproduced by the universal relation for any values 
of K. The disagreement of the theoretical prediction with the yield data is expected to 
be due to an unsuitable estimate of Uq, as is the same as C. Indeed, Roth, Bohdansky, & 
Ottenberger (1979) have suggested Uq = 8.5 eV for Al 2 03 from the behavior of sputtering 
yield at impact energies below the maximum yield by He + . Adopting this surface binding 
energy, we can get a better fit with K = 0.08 (Fig. le), although the agreement of the 
universal relation with the yield data is the worst among grain materials considered in this 
paper. 

To our knowledge, little or no experimental data on sputtering yield exist for MgO, FeS 
and Fe 3 04. Thus, we employ the results by the EDDY calculations to extract the values 
of K for these materials. In the figures from If to lh (If for MgO, lg for FeS, and lh for 
Fe 3 4 ), we show the EDDY results and the best-fitting theoretical curves. The agreements 
are not always good for H + and Ar + at low energies, which is also true for Si0 2 , but the 
universal relation well reproduces the yield data of Si0 2 . Thus, we consider that the fits are 
reasonable and satisfactory for these materials. The values of K determined are 0.06, 0.18 
and 0.15 for MgO, FeS and Fe 3 4 , respectively. 

Table 4 summarizes the value of K for each grain species evaluated by fitting to the 
available experimental data and/or the results of the EDDY simulations. We confirm that 
the derived value of K can well reproduce the experimental data of compounds as well as 
single-element materials, by adopting appropriate average values for the atomic and mass 
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numbers as proposed by Tielens et al. (1994). By slightly modifying the function a^fxi) in 
the universal relation, we realized that the calculated sputtering yields well agree with the 
available experimental data than before. The modified universal relation derived here could 
be applicable to evaluate the sputtering yield for any combinations of targets and projectiles 
of astrophysical interest as long as the experimental data are available. 



6. Physics of dust grains in shock 

The collisional interactions between dust and a hot gas in shock efficiently transfers 
charge, momentum and energy between the two phases. As a result, dust grains acquire 
the electric charge, undergo the resistance to their motion through the gas, are eroded by 
sputtering, and are heated by collisions with gas. Also, heated dust grains cool the gas 
through their thermal emission. 

Electric charge on a grain plays an important role in grain physics in hot plasmas 
through the Coulomb interaction and the Lorentz force in the presence of magnetic field. 
Here, we concentrate on the destruction of dust grains in fast (> 100 km s _1 ) and hot 
(> 10 5 K) nonradiative shocks. For T > 10 6 K, the effect of the electric charge of dust 
grain can be neglected because the dimensionless potential parameter <p defined by the ratio 
of the electric potential acquired by dust to the kinetic energy of gas is approximated by 
4> ~ 10 5 /T and is significantly smaller than unity (Draine & Salpeter 1979; McKee et al. 
1987). Also, the erosion rate of refractory grains by thermal sputtering quickly decreases for 
T < 10 6 K (Draine & Salpeter 1979). Correspondingly, the erosion rate of dust grans by 
non-thermal sputtering decreases quickly for the relative velocity of dust to gas < 200 
km s -1 . Therefore, it can be expected that the electric charge does not significantly affect 
the destruction of dust grains in nonradiative shocks considered in this paper. Also, it is 
assumed that the early ISM is not pervaded by a magnetic field because the magnetic field in 
the early universe is considered to be very weak (Gnedine, Ferrara, & Zweibel 2000). Thus, 
we formulate the basic equations describing the motion of dust grains and their erosion by 
sputtering in postshock flow, neglecting the effect of charge on dust grain. 



6.1. Dynamics of dust grain 

Once dust grains being at rest in the preshock ISM encounter the interstellar shock, 
they move ballistically behind the shock front with the initial velocity relative to gas w d ~ 
(3/4)V^hock where V^ock is the shock velocity, and are eroded by non-thermal sputtering. 
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A grain streaming through the ionized gas is decelerated by a drag force due to the direct 
collisions with gas particles, and thus the velocity of dust relative to gas decreases. The 
drag force acting on a grain results from the momentum transfer from gas to dust. Under 
the assumption that a gas particle reflects specularly from grain surface and that the kinetic 
temperature of gas is the same for all gas species, the deceleration rate of a spherical grain 
of radius a with the velocity relative to gas Wd is given by (e.g., Draine & Salpeter 1979) 



where nn is the number density of hydrogen atom, k is the Boltzmann constant, T is the gas 
temperature, pd is the bulk density of grain and sf = miw\j1kT with the mass of gas 
species %. The summation is taken over all gas species, and Ai is the number abundance of 
gas species % relative to that of hydrogen atom. For the function Gi(s{) whose exact formula 
is given by Baines, Williams, & Asebiomo (1965), the following analytical approximation 
has been proposed by Draine & Salpeter (1979) 



with the accuracy within 1 % for < Sj < oo. In the calculations, we adopt this formula. It 
can be seen from Equation (19), grains with smaller size and/or higher relative velocity can 
be effectively decelerated by drag force and come to comove with gas in shock. 



The destruction of dust grains in high- velocity and high-temperature nonradiative shock 
is dominated by sputtering. Actually, we can neglect the thermal evaporation of dust, 
because the temperature of dust cannot be raised as high as their sublimation temperature 
(~ 1000-2000 K) for the preshock gas density considered in this paper (e.g., Dwek 1987). 
Furthermore, we neglect the shattering process and partial vaporization caused by grain- 
grain collisions, since the number density of dust particles is very small in postshock flow 
and collisions between dust grains are extremely rare without a magnetic field. 

The erosion rate of grain by sputtering is calculated by taking the angle-averaged sput- 
tering yield (Y^E^q = 2Y l °(E i ) (Draine & Salpeter 1979) and is given by (e.g., Dwek et al. 
1996) 
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6.2. Grain destruction by sputtering 
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where q = Ei/kT and m sp is the average mass of the sputtered atoms. Equation (21) can 
be reduced to the following equations for the two extreme cases of (Dwek & Arendt 1992; 
Tielens et al. 1994); one is for Sj 1 where a stationary grain suffers the thermal sputtering 
caused by collisions originating from thermal motion of gas, and the erosion rate is written 
as 



The other is in the limit of Sj — > oo where the non-thermal sputtering erodes a hypersonic 
grain with the rate given by 



Figure 2 shows the erosion rate of each dust species by sputtering in the gas with the 
primordial elemental composition of metallicity Z = 10~ 4 Z Q given in Table 1; Figure 2a 
for thermal sputtering versus the gas temperature T from Equation (22) and Figure 2b for 
non-thermal sputtering versus the velocity of dust relative to gas Wd from Equation (23). 
The erosion rate by thermal (non-thermal) sputtering steeply increases from T ~ 10 5 K 
(w d ~ 30 km s^ 1 ), reaches a peak at T = 4-20 xlO 7 K (w d = 500-1300 km s" 1 ), and then 
slowly decreases with increasing T (w^). This behavior of the erosion rate as a function of 
T and Wd reflects the dependence of sputtering yield on impact energy (see 5.2). Among 
dust species considered in this paper, C grain has the lowest erosion rate at T > 2 x 10 6 K 
(wd > 200 km s _1 ), which is about one order of magnitude lower than that of FeS grain with 
the highest rate at T > 10 7 K (wd > 400 km s _1 ). For thermal sputtering, the erosion rate 
da/dt for the other dust species is ~ 1.2 x 10~ 6 n H A*m yr _1 cm 3 within a factor of three at 
T > 2 x 10 6 K. For non-thermal sputtering at Wd < 200 km s" 1 , the destruction by He + is 
dominant, and the erosion rate of AI2O3 is the lowest. 



The collisions with gas can also heat dust grains in postshock flow. The collisional 
heating rate H(a,T,nu) of dust grain with radius a is presented by Dwek & Arendt (1992) 

as 




(22) 
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6.3. Grain heating 
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where r)i(a, Cj) is the fraction of kinetic energy of gas species % deposited into the dust grain. 
The values of r)i(a,€i) are calculated by comparing an effective grain thickness a cS (4a/3 
for a spherical grain) with the stopping range of the incident particle l s . When a e g > 
l s , the incident particle deposits almost all of the energy into dust grains, otherwise, the 
difference between a c g and l s is exploited to determine the fraction of deposited energy. 
The approximations for the stopping range of ion species considered here are given in Dwek 
(1987). For the electron stopping range, we adopt the approximation derived by Tabata, 
Ito, & Okabe (1972) and Iskef, Cunningham, & Watt (1983). This electron stopping range 
reduces the heating rate by ~ 30 % at 10 6 K < T < 10 9 K compared with that given by 
Dwek (1987). The values and detailed derivation of ^(a, E) will be presented elsewhere 
(Nozawa et al. in preparation). 



7. Results of calculation of dust destruction 

In this section, we show the results of numerical simulations for dust destruction in non- 
radiative shock, focusing on the results calculated for the ISM-dust specified by the unmixed 
grain model described in § 3. In 7.1, we demonstrate the time evolution of temperature and 
density of gas and the destruction of dust grains in the interstellar shock that is driven by 
the SN model of C20 with E 5i = 1 and M pr = 20 M and is propagating into the ISM with 
n H,o — 1 cm" 3 . In 7.2, we investigate the dependences of the destruction efficiency and the 
mass of gas swept up by shock on E 5 i, M pr and nu,o, and derive the approximation formula 
for the time-scale of destruction for each dust species in the early universe. The effects of 
the cooling processes of gas on the evolution of nonradiative shock and the destruction of 
dust are briefly discussed in 7.3. 



7.1. Dust destruction in interstellar shock for C20 model 

7.1.1. Time evolution of temperature and density of gas in shock 

Figure 3 shows the structures of density (upper panel) and temperature of gas (lower 
panel) in nonradiative shock at given times, and the solid lines in Figure 4 give the time 
evolution of the shock velocity Vshock and the gas temperature T s h oc k at the shock front. 
As shown in Figure 3, the gas density increases to 4 times that in the ISM and the gas 
temperature steeply raises up at the shock front which is indicated by the arrow for a given 
time step. With the initial shock velocity ~ 6000 km s _1 , the gas temperature at the shock 
front remains above 10 8 K until t ~ 200 years, and both V s h oc k and T s h oc k decrease with time. 
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At t > 5 x 10 4 years after the explosion, the density of gas in the inner postshock region is 
about one order of magnitude lower than that in the unshocked region, and the low-density 
hot bubble with temperature of several times 10 6 K is formed (Fig. 3). In this model, V^hock 
decelerates below 100 km s _1 at t tT ~ 10 5 years when the shock front travels the distance of 
~ 30 pc and T s h oc k drops down to several 10 4 K. After then, a dense thin-shell forms at the 
shock front because of the effective line cooling, and the SNR enters into radiative phase. 

7.1.2. Dynamics and destruction of dust grains in shock 

The motion of dust grains (upper panel) and the evolution of their sizes (lower panel) in 
postshock flow as a function of time are depicted in Figure 5, as an example for C grains with 
the initial size of 0.01 /j,m (dashed), 0.1 /im (dotted), and 1 /im (solid). The thick solid curve 
in the upper panel indicates the position of the shock front. As dust grains initially at rest in 
the ISM intrude into the blast wave, they depict the different trajectories depending on their 
initial sizes, which clearly demonstrates that the dust grains are segregated and subjected 
to the different sputtering processes in the postshock flow; small grains with initial radius 
of 0.01 /im are quickly decelerated by drag force, are trapped into gas near the shock front, 
and are completely destroyed by thermal sputtering. The 0.1 /im-sized grains are gradually 
decelerated to comove with gas in 7 x 10 3 years after entering into the shock. Large-sized 
grains with radius of 1 /im are not so much suffered from deceleration by drag force and 
continue to keep the high velocity relative to gas. As a result, they are subjected to non- 
thermal sputtering, but their sizes are reduced very little, partly because they stay in the 
inner region of postshock flow where the gas density is lower than that near the shock front 
(see Fig. 3). 

The modification of size distribution of each dust species due to destruction by sputtering 
is illustrated in Figure 6; Figure 6a for the initial size distribution before destruction and 
Figure 6b for the size distribution after destruction. Since the erosion rate by sputtering 
does not strongly depend on grain size (see 6.2), small grains are predominantly destroyed 
regardless of grain species; the number of small-sized grain such as A1 2 3 is greatly reduced. 
C, SiC>2 and Fe grains whose initial size distributions are lognormal with the relatively large 
average size are eroded but not completely destructed, and increase the numbers of smaller 
ones. Grains with the size larger than 0.1 /im are little affected by erosion. Note that the 
size distribution summed up over all grains species gets flatter for the smaller radius with 
time compared with the corresponding initial size distribution approximated by a power-law 
formula with index of —2.5, while that for the larger radius remains almost unchanged. 



7.1.3. Efficiency of dust destruction 



The destruction efficiency 6j of each grain species is given in Table 5 (see unmixed 
grain model) along with the initial average radius. It is expected that the efficiency of dust 
destruction are higher for grain species with the smaller average size because the smaller 
grains are predominantly destructed by sputtering. In fact, the efficiency of destruction of 
A1 2 3 grain with the smallest average radius is 0.667 and is the highest among dust species 
considered here. Si and Fe grains with initial average radii larger than 0.1 /xm have the 
destruction efficiencies less than 0.2; 0.13 for Si grain and 0.15 for Fe grain. However, the 
destruction efficiency of C grain with the smaller average radius is smaller than that of Si0 2 
grain because of the lowest erosion rate among all grain species at T > 2 x 10 6 K (see Fig. 
2). Also, FeS grain has higher destruction efficiency (0.578) than those of Mg2SiC>4 (0.451) 
and MgO grains (0.505) despite of the larger initial average radius, which reflects not only 
the higher erosion rate at T > 10 7 K but also the mass distribution much more weighted 
toward the smaller grains. Likewise, the destruction efficiency of MgSi0 3 grain with average 
size comparable to that of Mg 2 Si0 4 and MgO is high (0.637) because of the lack of large 
grains. These facts indicate that the efficiency of dust destruction depends on the initial size 
distribution of dust grains as well as on the sputtering yield. 

In order to clarify the effect of the initial size distribution on the efficiency of dust 
destruction, we present the results of dust destruction for the mixed grain model calculated 
with the same SN model and value of %,o as those for the unmixed grain model. Figure 7 
shows the size distribution of each grain species (AI2O3, MgSiC>3, Mg2Si04, SiC>2 and FeaC^) 
before destruction (Fig. 7a) and after destruction (Fig. 7b), and the destruction efficiencies of 
these grains are tabulated in Table 5 (see mixed grain model). As is the same as the unmixed 
grain model, the numbers of A1 2 3 and Fe 3 4 grains with the small average size (< a few 
tens A) are considerably reduced, which results in high efficiencies of destruction (> 0.7). 
For MgSi0 3 , Mg 2 Si0 4 and Si0 2 grains with lognormal size distributions, the erosion of large 
grains leads to the increase of the number of smaller ones, and their destruction efficiencies 
span the range of 0.4 < tj < 0.64. Note that although the average size is twice larger than 
that in the unmixed case, the destruction efficiency (0.59) of Mg2Si04 grain in the mixed 
case is significantly larger than that (0.45) in the unmixed case. The reason is that Mg2SiC>4 
grains formed in the unmixed ejecta have power-law-like size distribution and include the 
large-sized grains of > 0.1 almost undestructed in shock; the average size is not always 
suitable for assessing the feasibility of dust destruction. Thus, we conclude that the dust 
destruction efficiency is very sensitive to the initial size distribution. The mass fraction of 
dust destructed reaches up to 34 % for the unmixed grain model and 50 % for the mixed 
grain model. 
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7.2. Time-scale of dust destruction in the early universe 

In this subsection, we investigate the dependences of the efficiency of dust destruction €j 
and the mass of gas swept up by shock M swep t on £51, M pr and nu,o, and derive an analytic 
formula describing the time-scale of dust destruction for the unmixed grain model. Figure 8 
shows the destruction efficiency of each grain species versus SN explosion energy calculated 
for n H ,o = 1 cm -3 ; Figure 8a is for A1 2 3 , FeS, Mg 2 Si0 4 and Fe grains, and Figure 8b is for 
MgSi0 3 , MgO, Si0 2 , C and Si grains. Also, the overall efficiency of dust destruction, which 
is defined as the ratio of the total mass of dust destructed to the total mass of dust swept up 
by shock, is plotted in Figure 8a. The SN models used for the calculations are distinguished 
by open circles (CCSNe), open squares (HNe) and filled triangles (PISNe). 

The destruction efficiencies for each grain species are almost the same for CCSNe and 
a PISN with the explosion energy of E 5 i = 1 irrespective of the progenitor mass M pr , and 
increase with increasing E 5 i. Note that high explosion energy with E 5 i > 10 causes the 
temperature of gas in postshock flow to raise up as high as 10 9 K, but this does not directly 
influence on the efficiency of dust destruction because the increase of gas temperature does 
not always lead to the enhancement of the erosion rate by sputtering (see Fig. 2a). The 
reason for the increased efficiency with increasing E 5 i is considered as follows. The high- 
velocity shock (the initial shock velocity > 10 4 km s _1 ) generated by the energetic SN 
explosion induces high velocity of dust relative to gas. Then, dust grains are efficiently 
decelerated by drag force (see Equation (19)), are trapped in the high-density region near 
the shock front, and are significantly eroded by thermal sputtering. Furthermore, since the 
high shock velocity takes much longer time to drop down to 100 km s _1 (see Fig. 9), dust 
grains are immersed in a hot plasma for a long time, which causes even larger grains to be 
more eroded. 

It should be pointed out here that the efficiencies of destruction for the models of PI 70 
with E 51 = 20 and P200 with E 51 = 28 are a little higher than that for H30 with E 51 = 30. 
The reason is considered that the ejecta mass of PISNe is more than 6 times larger than 
that of HNe, and the resulting longer duration of the free expansion phase causes the longer 
truncation time than HNe (Fig. 9). On the other hand, the efficiency of destruction for 
P150 with M pr = 150 M and E 5 i = 1 is almost the same as that for CCSNe with E 5 i = 1 
because the initial shock velocity (~ 3000 km s _1 ) much lower than that of CCSNe makes 
the truncation time comparable to CCSNe. Therefore, although the destruction efficiency 
6j is almost independent of M pr as long as £51 = 1, tj for PISNe higher than that for HNe 
at £51 > 10 reflects the difference in the explosion mechanism depending on the progenitor 
mass. 

To examine the dependence of ej on E 51 for each type of SN, we calculate the coefficients 
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a hj and b ld in Equation (Al) given in Appendix for SNe II (C13, C20, C25, C30, H25 and 
H30) and PISNe (P150, P170 and P200) separately, and the calculated overall efficiencies 
of dust destruction are indicated by dotted (SNe II) and dashed lines (PISNe) in Figure 
8a. Although the difference in the efficiency of dust destruction between SNe II and PISNe 
increases with increasing E51, the deviations from the values calculated by a±j and bij for 
all SN models are at most about 10 % at £51 = 30. Thus, we consider that ej is almost 
independent of M pr as long as E 51 is the same. For reference, we tabulate the values of a 1: j 
and b lt j for SNe II and PISNe in Table 6. The mass of gas swept up by shock M swept and 
the truncation time t tT calculated for nn,o = 1 cm" 3 are presented in Figure 9 as a function 
of SN explosion energy. As is the same as € j: M swcpt and t tI are almost the same for E 51 = 1 
regardless of M pr , and increase with increasing E 5 i. 

Next, we show the dependences of €j and M swcpt on the preshock gas density 71h,o- 
Figure 10a plots the overall efficiency of dust destruction versus nn,o for C20 (crosses), H25 
(open circles) and H30 models (filled triangles) with £51 = 1, 10 and 30, respectively. The 
efficiency of dust destruction increases with increasing nn,o, since higher gas density causes 
more frequently collisions between dust and gas to erode efficiently the surface of dust grains 
by sputtering; for example, the mass fraction of dust destroyed for H30 model is 78 % for 
nn,o = 10 cm -3 , but only 23 % for nn,o = 0.1 cm~ 3 . In Figure 10b, we present the mass 
of gas swept up by shock M swep t as a function of nn,o for C20, H25 and H30 models. Note 
that M swept decreases with increasing nn,o, because shock wave more quickly decelerates and 
travels only small distance. 

The approximation formulae presented in Appendix being combined, the time-scale of 
destruction for each dust species by the interstellar shock in the early universe is presented 
by 

-1 (T? n 4144£/ 51 n H,0 M & , n r\ 

T SN,i = e i(^5i, n H ,o) T7 7sn, (25) 

for all SN models, where the dependences of tj on E 51 and n H ,o are given by Equations 
(Al) and (A3), respectively. We derived the time-scale of dust destruction as a function of 
not only the explosion energy of SN but also the gas density in the ISM. The swept-up gas 
mass is proportional not to £51 but to ~ E®f, being different from the formula proposed 
by McKee (1989). The time-scale of dust destruction derived here could be applicable to 
investigate the time evolution of dust grains in the early universe. Especially, the difference 
of the efficiency of dust destruction for each grain species may have great influences on the 
amount and size distribution of dust grains residing in the early interstellar space. 
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7.3. The effects of the cooling processes of gas on dust destruction 

Figure 1 1 shows the cumulative energy lost by the atomic process (dashed) , the inverse 
Compton scattering (dotted) and the thermal emission from dust (solid) calculated for the 
SN model of C20 and for the ISM with parameters of n H ,o = 1 and Z = 10~ 4 Z Q at the 
redshift of z = 20. The total energy lost by these cooling processes is drawn by the thick 
solid curve. The inverse Compton cooling is comparable to that by the atomic process in 
the early phase of the SNR (t < 10 4 years). As the gas temperature decreases, the atomic 
process by H and He line coolings becomes dominant. Compared with the above two cooling 
processes, the thermal emission from dust is extremely low and contributes only less than 0.1 
% to the total energy loss. The transition of nonradiative shock to radiative shock occurs 
when the cumulative energy loss reaches to 0.01 % of the explosion energy. Thus, only 
the H and He line coolings affect the time evolution of nonradiative shock and the dust 
destruction efficiencies in the ISM with metallicity less than Z = 10~ 4 Z Q corresponding to 
the dust-to-gas mass ratio of 4.5 x 10~ 7 in this paper. 

However, it is expected that the effects of cooling of gas by dust come to be important 
for the ISM with higher dust-to-gas mass ratio. The result of the calculation for dust 
destruction in the ISM with Z = 10 -4 Z Q and the dust-to-gas mass ratio 4.5 x 10 -3 shows 
that the truncation time and the overall efficiency of dust destruction decrease by 4.4 % and 
2.6 %, respectively, compared with the results for dust-to-gas mass ratio of 4.5 x 10~ 7 , though 
the gas temperature in postshock gas decreases by about 20%. Even if the dust-to-gas mass 
ratio is raised up to by four orders of magnitude, the cooling of gas by dust in nonradiative 
shock does not have the significant effects on the efficiency of dust destruction as well as the 
evolution of nonradiative shock. 

Furthermore, the inverse Compton cooling has influences on the evolution of nonradia- 
tive shock in the early universe. Because the cooling rate by the inverse Compton scattering 
is proportional to n e T(l + z) A and that by atomic process is proportional to n e nn ~ n^, the 
contribution of the inverse Compton cooling is enhanced in postshock flow with the lower gas 
density and/or the higher gas temperature at the higher redshift. At the redshift of z — 40, 
for the SN model C20 with the preshock gas density of nn,o = 1 cm -3 , the truncation time 
is only 2 % shorter than that calculated at z — 20. Thus, in this case, the inverse Compton 
cooling seems not to affect the destruction efficiency of dust grains. However, more system- 
atic studies covering wide ranges of the SN explosion energy and the gas density in the ISM 
are necessary to reveal the effect of the inverse Compton cooling. 

Note that the cooling function of gas for Z = 10~ 3 Z Q is almost the same as that for the 
zero metal case (Sutherland & Dopita 1993). Thus, the parameters of the ISM considered 
in this paper (Z = 10~ 4 Z Q and z = 20) is competent for investigating dust destruction in 



-25 - 



the early universe, and the time-scale of dust destruction derived in 7.2 is applicable over 
the wide ranges of metallicity of gas in the ISM (Z < 10~ 3 Z & ) and redshift (z < 40). 

8. Summary 

We investigate the destruction of dust grains in the interstellar shocks driven by SNe as 
the second step to reveal the evolution of dust in the early universe, based on the dust models 
obtained by Nozawa et al. (2003). We focus on dust destruction in nonradiative shock where 
dust grains are predominantly destroyed by non-thermal and thermal sputtering because of 
high temperature and high velocity of gas. The sputtering yields for the combinations of 
dust and ion species of interest to us are evaluated by applying the universal relation for 
sputtering yield with a slight modification and by determining the value of K by fitting to 
the available experimental data and/or the results of the EDDY simulations. The modified 
universal relation derived here can present the better fits to the available experimental data 
than the previous one. 

In the calculations of dust destruction, we solve the time evolution of gas temperature 
and density in spherically symmetric shocks, adopting the hydrodynamical models by Umeda 
& Nomoto (2002) as the initial condition for interstellar shock and including the cooling of 
gas by the atomic process, the inverse Compton scattering and the thermal emission of dust. 
The erosion of dust by thermal and non-thermal sputtering caused by motion of dust relative 
to gas is carefully treated by taking into account the size distribution of each dust species. 
The results of calculations are summarized as follows. 

(1) Because the sputtering predominantly destroys the small grain, the number of small- 
sized grain such as A1 2 3 and MgSi0 3 is greatly reduced. The erosion of C, Si0 2 and Fe 
grains whose size distributions are lognormal-like with the average size larger than 0.01 /iin 
increases the numbers of smaller ones. The size distribution summed up over all grains 
species becomes flatter for the small radius compared with the initial size distribution, while 
that for the radius larger than 0.2 /im remains almost unchanged. 

(2) The efficiency of dust destruction is higher for the grains with the small average size 
such as AI2O3 and MgSiOa and with the power-law-like size distribution such as FeS, Mg 2 Si04 
and MgO. On the other hand, Si, Fe, C and SiC>2 grains, which have the lognormal-like size 
distribution and the relatively large average radius, have lower efficiency of destruction. The 
detailed analysis of the behavior of dust destruction efficiency for each grain species indicates 
that not only sputtering yields but also the initial size distribution plays a crutial role in the 
efficiency of dust destruction by sputtering. 
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(3) The efficiency of destruction €j for each dust species increases with increasing the 
explosion energy £ 51 and/or with increasing the preshock gas density »h,0) but is almost 
independent of the progenitor mass M pr of SN as long as £51 is the same. The destruction 
efficiency €j as a function of £51 is well reproduced by a power-law formula given by Equation 
(Al). The dependence of tj on nn,o is well expressed by the quadratic equation (A3) in terms 
of log(n H ,o)- 

(4) As is the same as 6j, the mass of gas swept up by shock wave M swept is the increasing 
function of £51 and is approximated by a power- law formula given by Equation (A2). How- 
ever, M swep t decreases with increasing the gas density in the ISM, and its dependence is also 
reproduced by a power-law formula whose index is given by Equation (A4) as a function of 
£ 51 . Finally, by combining these results, we present the analytic formula for the time-scale 
of destruction for each grain species in the early universe as a function of £51 and tih,o- The 
derived time-scale of dust destruction can be employed to investigate the time evolution of 
the amount of dust grains in the early universe. 

(5) In the early universe, only the H and He line coolings affect the time evolution of 
nonradiative shock and the efficiencies of dust destruction. The thermal emission from dust 
grains is not so much important even for the ISM with the dust-to-gas mass ratio of ~ 1CT 3 
as long as the metallicity in the ISM is Z < 10~ 3 Z Q . Also, it is expected that the inverse 
Compton cooling does not have the great effects on the evolution of shock at high redshifts. 
Thus, the time-scale of dust destruction derived in this paper is applicable for ranges of 
metallicity of gas in the ISM of Z < 10~ 3 Z Q and redshift of z < 40. 

In this study, we focused on dust destruction in nonradiative shock where dust grains in 
the ISM are considered to be predominantly destroyed. We should be mentioned here that a 
part of dust grains formed in the ejecta of SNe are expected to be destroyed by reverse shock 
penetrating into the ejecta. Because sputtering is the most dominant destruction process of 
dust grains in reverse shock as well, the simulation code of dust destruction constructed in 
this paper can be applied to explore the destruction of dust grains by reverse shocks. This 
work is now in progress. Furthermore, in radiative shock with the shock velocities below 100 
km s -1 , shattering and partial vaporization by grain-grain collisions become important in 
the presence of a magnetic field. A magnetic field frozen into the gas causes charged grains 
to gyrate and accelerate behind the shock front and results in the enhancement of not only 
the sputtering rate but also grain-grain collision frequency. The shattering by grain-grain 
collision mainly leads to the redistribution of grain size (Tielens et al. 1994; Jones, Tielens, 
& Hollenbach 1996). The sophisticated study on this subject is left for the future work. 

The authors are grateful to the anonymous referee for critical comments which have 
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part by a Grant-in-Aid for Scientific Research from the Japan Society for the Promotion of 
Sciences (13640229). 



A. The appriximation formulae for ej and M swept 

In this Appendix, we present the approximate formulae for the efficiency of dust de- 
struction 6j and the mass of gas swept up by shock M swcpt as function of the SN explosion 
energy E 51 and the gas density n H ,o i n the ISM. The energy dependence of the destruction 
efficiency for each grain species can be well reproduced by a power-law formula given by 

€j = a 1;j E b 5 \'\ (Al) 

The solid lines in Figure 8 indicate the least-squares fittings by this power-law formula for 
all SNe whose numerical coefficients aij and bij are tabulated in Table 6. The dependence 
of M swcpt on E 5 i is well fitted by a power-law formula 

M swept = a 2 El\ (A2) 

using the coefficients a 2 and b 2 given in Table 6, and is depicted in Figure 9 for all SNe 
(solid), SNe II (dotted), and PISNe (dashed). The destruction efficiency ej for each grains 
species j as well as the overall destruction efficiency as a function of n H) o is expressed by 

loge,- = Cj (logn H ,o) 2 + dj logn H ,o + e j? (A3) 

and the coefficients Cj, dj and tj derived by the least squares fittings are given in Table 7 for 
the models C20, H25 and H30, respectively, along with the coefficients for overall destruction 
efficiency. The solid curves in Figure 10a depict the overall destruction efficiencies calculated 
by Equation (A3) for the models C20, H25 and H30. This formula reproduces the destruction 
efficiencies calculated by the simulations with the accuracy less than 10 % for n H ,o < 10 cur 3 . 
The calculated M swept can be reproduced by a power-law formula M swept cx nf^o where the 
index g is weakly dependent on E 5 i and is approximated by 

g = -0.U2E™ 63 . (A4) 

This formula gives good agreements with the results of simulations, as represented by the 
solid lines in Figure 10b. 
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Table 1. Elemental Abundances of Gas in the ISM 



z/z Q 




Ac 


^4n 


Ao 


A e 


1(T 4 


8.04 x 1(T 2 


3.62 x 10~ 8 


1.12 x 10~ 8 


2.70 x 10~ 7 


1.128 



Note. — The number abundance A^ relative to that of hydrogen atom 
in the ISM. Note that the abundances of C, N and O are calculated by 
multiplying the metallicity of Z — 10~ 4 Z Q and the primordial ratios 
in Table 4 of Sutherland & Dopita (1993). 
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Table 2. Dust Models Used in the Calculations 







unmixed j 


jrain model 


mixed 


grain model 


dust species 


Pd,j (g cm 3 ) 


Aj a 


avc,i (A*m) 


Aj 


a ave ,j (/Jin) 


C 


2.26 


0.103 


0.029 






Si 


2.32 


0.304 


0.250 






Fe 


7.89 


0.050 


0.148 






FeS 


4.84 


0.156 


0.0088 






A1 2 3 


3.98 


0.002 


0.0007 


0.001 


0.0003 


MgSiOg 


3.18 


0.010 


0.0045 


0.102 


0.011 


Mg 2 Si0 4 


3.20 


0.219 


0.0066 


0.205 


0.010 


Si0 2 


2.64 


0.100 


0.047 


0.588 


0.033 


MgO 


3.56 


0.056 


0.0058 






Fe 3 4 


5.21 






0.103 


0.0022 



Note. - The models of dust grains are adopted from the results of the 
calculations of dust formation in the ejecta of Pop III SNe with M pr = 20 
M & by Nozawa et al. (2003). The bulk density of grain species j is denoted 
by • The mass fraction Aj and the average size a avc j of each grain species 
are given for the unmixed and mixed grain models. Note that the average 
radius of Si grain in the unmixed grain model is calculated for the grains 
larger than 0.03 fxm which occupy 99.8 % of the total mass. 
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Table 3. Models of Supernova Ejecta 





Progenitor mass 


Explosion energy 


Model 


(M e ) 


E 51 (10 51 ergs) 


C13 


13 


1 


C20 


20 


1 


C25 


25 


1 


C30 


30 


1 


H25 


25 


10 


H30 


30 


30 


P150 


150 


1 


P170 


170 


20 


P200 


200 


28 



Note. - The models of the SN ejecta are 
adopted from Umeda & Nomoto (2002). In the 
model names, the labels C, H and P represent 
CCSNe, HNe and PISNe, respectively, and the 
numerical value denotes the mass of the progeni- 
tor in units of solar mass. 
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Table 4. Parameters for Calculations of Sputtering Yield and References for Yield Data 



dust species 


Uq (eV) 


z d 


M d 


K 


references 


C 


4.0 


6 


12 


0.61 


1, 2, 3, 4, 5 


Si 


4.66 


14 


28 


0.43 


3, 4, 6, 7, 8, 9, 10, 11 


Fe 


4.31 


26 


56 


0.23 


3, 4, 5, 7, 8, 12, 13 


FeS 


4.12 


21 


44 


0.18 


EDDY 


A1 2 3 


8.5 


10 


20.4 


0.08 


3, 14, 15 


MgSiOg 


6.0 


10 


20 


0.1 


16 


Mg 2 Si0 4 


5.8 


10 


20 


0.1 


16 


Si0 2 


6.42 


10 


20 


0.1 


3, 14, 15, 17, 18 


MgO 


5.17 


10 


20 


0.06 


14, 15, EDDY 


Fe 3 4 


4.98 


15.7 


33.1 


0.15 


EDDY 



Note. — The surface binding energy is represented by Uo in units of 
eV. The mean atomic number and mean mass number of the target 
material are and Md, respectively, and the free parameter K is 
determined by fitting to the available experimental data and/or the 
results of the EDDY simulation. See text for details. 

References. — (1) Bohdansky, Bay, & Ottenberger (1978); (2) Bo- 
hdansky, Roth, & Sinha (1976); (3) Roth, Bohdansky, & Ottenberger 
(1979); (4) Rosenberg & Wehner (1962); (5) Hechtl, Bohdansky, & 
Roth (1981); (6) Laegreid & Wehner (1961); (7) Southern, Willis, & 
Robinson (1963); (8) Eernisse (1971); (9) Sommerfeldt, Mashkova, 
& Molchanov (1972); (10) Coburn, Winters, & Chuang (1977); (11) 
Blank & Wittmaack (1979); (12) von Seefeld et al. (1976); (13) Bo- 
hdansky, Bay, & Roth (1977); (14) Bach (1970); (15) Nenadovic et al. 
(1990); (16) Tielens et al. (1994); (17) Bach, Kitzmann, & Schroder 
(1974); (18) Edwin (1973) 
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Table 5. The Efficiency of Dust Destruction for C20 Model with n H ,o = 1 





unmixed 


grain model 
a avCii (/im) 


mixed 

e i 


grain model 
a ave ,j (/im) 


c 


0.247 


0.029 






Si 


0.134 


0.250 






Fe 


0.154 


0.148 






FeS 


0.578 


0.0088 






A1 2 3 


0.667 


0.0007 


0.794 


0.0003 


MgSi0 3 


0.637 


0.0045 


0.631 


0.011 


Mg 2 Si0 4 


0.451 


0.0066 


0.586 


0.010 


Si0 2 


0.411 


0.047 


0.399 


0.033 


MgO 


0.505 


0.0058 






Fe 3 4 






0.741 


0.0022 


overall 


0.340 




0.497 
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Table 6. Coefficients for the Approximation Formula (Al) and (A2) 





for all SNe 




for SNe II 




for PISNe 


dust species 














a i,j 




C 


2.532E-01 


1.660E- 


■01 


2.526E-01 


1.407E- 


■01 


2.580E-01 


1.820E-01 


Si 


1.380E-01 


2.371E- 


■01 


1.377E-01 


1.979E- 


■01 


1.414E-01 


2.631E-01 


Fe 


1.583E-01 


2.402E- 


■01 


1.586E-01 


2.005E- 


■01 


1.594E-01 


2.707E-01 


FeS 


5.842E-01 


7.897E- 


•02 


5.850E-01 


6.758E- 


•02 


5.834E-01 


8.857E-02 


A1 2 3 


6.724E-01 


5.649E- 


■02 


6.733E-01 


4.803E- 


•02 


6.709E-01 


6.393E-02 


MgSiOs 


6.429E-01 


6.358E- 


•02 


6.438E-01 


5.421E- 


•02 


6.417E-01 


7.165E-02 


Mg 2 Si0 4 


4.585E-01 


1.146E- 


•01 


4.590E-01 


9.827E- 


•02 


4.595E-01 


1.275E-01 


Si0 2 


4.183E-01 


1.238E- 


■01 


4.189E-01 


1.054E- 


■01 


4.19E-01 


1.383E-01 


MgO 


5.117E-01 


1.046E- 


■01 


5.122E-01 


9.020E- 


•02 


5.124E-01 


1.160E-01 


overall 


3.458E-01 


1.307E- 


■01 


3.460E-01 


1.099E- 


■01 


3.476E-01 


1.463E-01 






a 2 


b 2 




a 2 


b 2 




a 2 


b 2 



-^swept 



4.144E+03 7.967E-01 4.118E+03 6.627E-01 4.481E+03 8.864E-01 
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Table 7. Coefficients for the Approximation (A3) 







C20 






H25 






H30 




dust species 




dj 


e i 




dj 






dj 


e i 


C 


-0.062 


0.41 


-0.61 


-0.07 


0.36 


-0.46 


-0.08 


0.33 


-0.39 


Si 


-0.044 


0.54 


-0.87 


-0.06 


0.50 


-0.67 


-0.072 


0.47 


-0.58 


Fe 


-0.057 


0.55 


-0.81 


-0.078 


0.49 


-0.61 


-0.096 


0.46 


-0.51 


FeS 


-0.088 


0.23 


-0.23 


-0.077 


0.18 


-0.16 


-0.073 


0.16 


-0.13 


A1 2 3 


-0.031 


0.14 


-0.17 


-0.030 


0.12 


-0.13 


-0.032 


0.11 


-0.11 


MgSiOs 


-0.060 


0.18 


-0.19 


-0.057 


0.14 


-0.13 


-0.055 


0.12 


-0.11 


Mg 2 Si0 4 


-0.087 


0.31 


-0.35 


-0.092 


0.25 


-0.24 


-0.097 


0.22 


-0.19 


Si0 2 


-0.086 


0.33 


-0.39 


-0.093 


0.27 


-0.27 


-0.097 


0.24 


-0.22 


MgO 


-0.079 


0.29 


-0.30 


-0.086 


0.23 


-0.20 


-0.089 


0.20 


-0.16 


overall 


-0.067 


0.33 


-0.47 


-0.067 


0.29 


-0.35 


-0.072 


0.26 


-0.30 
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Fig. 1. — Sputtering yields of each grain species versus incident energy of projectiles, (a) C, (b) 
Si, (c) Fe, (d) Si0 2 , (e) A1 2 3 , (f) MgO, (g) FeS, and (h) Fe 3 4 . The incident ion species are H+, 
D + , He + , Ne + and Ar + . The experimental data on sputtering yield are represented by open circles 
(H + ), asterisks (D + ), open squares (He + ), crosses (Ne + ) and open triangles (Ar + ), and the results 
of sputtering yield calculated by the EDDY code are denoted by filled circles (H + ), filled squares 
(He + ) and filled triangles (Ar + ). The solid curves show the best-fitting theoretical yields calculated 
by the universal relation. In (b) and (c), the dotted curves show the theoretical sputtering yields 
by Ar + projectile calculated by adopting cti given by Equation (17). See the electronic edition of 
the Journal for a color version of this figure. 
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Fig. 2. — The erosion rate of each dust species by sputtering in units of /im yr" 1 cm 3 
calculated for the elemental composition of gas with Z = 10~ 4 Z Q given in Table 1; (a) 
the thermal sputtering calculated by Equation (22) as a function of gas temperature T and 
(b) the non-thermal sputtering calculated by Equation (23) as a function of velocity of dust 
relative to gas w&. See the electronic edition of the Journal for a color version of this figure. 
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Fig. 3. — The structures of density (upper panel) and temperature of gas (lower panel) at 
given times in the interstellar shock that is driven by the SN model of C20 with Em = 1 and 
M pr = 20 Mq and is propagating into the ISM with nn,o — 1 cm -3 . The arrow in the lower 
panel indicates the position of the shock front at a given time. See the electronic edition of 
the Journal for a color version of this figure. 
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Fig. 4. — The time evolution of the shock velocity V^hock and the gas temperature T shock at 
the shock front for model C20 with tih,o = 1 cm -3 . See the electronic edition of the Journal 
for a color version of this figure. 
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Fig. 5. — The time evolution of position (upper panel) and size (lower panel) of dust grain 
in postshock flow for model C20 with n^o — 1 cm -3 . Dust species considered is C grain 
with the size of 0.01 jum (dashed), 0.1 /zm (dotted) and 1 /im (solid). The thick solid curve 
in the upper panel denotes the position of the shock front -R s hock- See the electronic edition 
of the Journal for a color version of this figure. 
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Fig. 6. — The size distribution of each dust species for the unmixed grain model; (a) for the 
initial size distribution before destruction and (b) for the size distribution obtained by the 
calculation of dust destruction for model C20 with ?t,h,o = 1 cm -3 . See the electronic edition 
of the Journal for a color version of this figure. 
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Fig. 7. — The size distribution of each dust species for the mixed grain model; (a) for the 
initial size distribution before destruction and (b) for the size distribution obtained by the 
calculation of dust destruction for model C20 with riH,o = 1 cm -3 . See the electronic edition 
of the Journal for a color version of this figure. 
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Fig. 8. — The efficiency of destruction of each grain species in the unmixed grain model calculated 
for riH,o = 1 cm -3 as a function of SN explosion energy; (a) for AI2O3, FeS, Mg2SiC>4 and Fe grains 
and (b) for MgSiC>3, MgO, Si02, C and Si grains. The SN models used for the calcularions are 
distinguished by open circles (CCSNe), open squares (HNe) and filled triangles (PISNe). Also, in 
(a), the overall efficiency of dust destruction is plotted. The linear solid lines for each grains species 
are the results calculated by the power-law formula given by Equation (Al) for all SNe. The dotted 
and dashed lines for the overall efficiencies of dust destruction are the results of calculations for 
SNe II and for PISNe, respectively. See the electronic edition of the Journal for a color version of 
this figure. 
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Fig. 9.— The mass of gas swept up by shock M swe pt and the truncation time ttr versus 
the explosion energy of SN for nH,o = 1 cm" 3 . The results for CCSNe, HNe and PISNe 
are represented by open circles, open squares and filled triangles, respectively. The linear 
lines are the power-law formula approximated by Equation (A2) for all SNe (solid), SNe II 
(dotted) and PISNe (dashed). See the electronic edition of the Journal for a color version 
of this figure. 
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Fig. 10. — (a) The overall efficiency of dust destruction versus tzh,o for C20 (crosses), H25 
(open circles) and H30 models (filled triangles). The solid curves are calculated by the 
approximation formula by Equation (A3), (b) The mass of gas swept up by shock for C20 
(crosses), H25 (open circles) and H30 models (filled triangles) as a function of riH,o- The linear 
lines are the power-law approximation formula of M swept oc n^ with g = —0.14:2E®f 63 . See 
the electronic edition of the Journal for a color version of this figure. 
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Fig. 11. — The time evolution of cumulative energy lost by the atomic process (dashed), the 
inverse Compton cooling (dotted) and the thermal emission from all dust grains (solid) in 
units of E$i for the SN model C20, riH,o — 1 an d Z = 10 -4 Z Q at the redshift of z = 20. 
The total energy loss is indicated by the thick solid curve. See the electronic edition of the 
Journal for a color version of this figure. 



